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We analyze the model of an atomic Bose-Einstein condensate (BEC) parametrically coupled to a 
molecular BEC via a photoassociation process. We show that an interplay of nonlinear inter- and 
intra-species interactions leads to the formation of mutually trapped states of a hybrid condensate, 
which are spatially localized even without a trap. The untrapped atomic-molecular condensates with 
a large number of atoms in either fraction are shown to be dynamically stable in the framework of 
a mean field theory. 
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PACS numbers: 03.75.Fi, 03.75.-b, 03.65.Ge, 33.80 Ps 



Recent progress in the field of the Bose-Einstein con- 
densation is associated with the possibility of creat- 
ing multi-component condensates of different, coher- 
ently coupled, atomic species. Furthermore, the pos- 
sibility of creating a hybrid condensate of atomic and 
molecular fractions has been suggested (l), and re- 
cently received some experimental support In 
such an atomic-molecular condensate, atoms would 
bind into molecules coherently, through a (reversible) 
photoassociation-induced stimulated recombination in 
two-body collisions. 

The problem of generation of coupled atomic and 
molecular condensates is fascinating from different points 
of view (see, e.g., Refs. [|f||)- Apart from the prospect 
of observing the first molecular condensate, this prob- 
lem offers a unique atom-optics analogue to the pro- 
cess of second- harmonic generation in nonlinear optics. 
The nonlincarity produced by coherent atomic- molecular 
coupling is expected to be responsible for a number of 
dynamical effects, including inter-species Josephson-like 
population oscillations ^|J|J^] , and "clumping" of the con- 
densate due to the atom-optical analogue of modula- 
tional instability In a yet to be realized atomic- 

molecular BEC (AMBEC), the atom- molecular para- 
metric coupling would compete with the atom-atom, 
molecule-molecule, and atom-molecule collisional inter- 
actions, which corresponds to the competing quadratic 
and cubic intensity-dependent response in nonlinear op- 
tics H , and suggests the existence of liquid- like droplets 
of the atomic-molecular condensates || . 

In this Letter, we present a comprehensive analysis of 
the mean-field model for an AMBEC created via a coher- 
ent stimulated Raman photoassociation process. Using 
the recently reported data for the intra- and inter-species 
scattering lengths 0, we find the existence domain for 
the stationary states of an AMBEC in an experimentally 
relevant parameter region. We then show that, under 
certain conditions, stationary states of an AMBEC are 
dynamically stable, which enables the creation of self- 
confined two-species three-dimensional AMBEC droplets, 
after the condensate is released from a confining trap. 
Untrapped droplets of the coherent matter waves resem- 



ble self-trapped spatial solitons supported by parametric 
interaction of optical waves || . 

We study the dynamics of the atomic and molec- 
ular condensates in a spherically symmetric, three- 
dimensional trap, dynamically coupled by a free-bound 
Raman photoassociation process in the framework 

of the coupled mean-field Gross-Pitaevskii (GP) equa- 
tions for the macroscopic wave functions of the two 
species, 
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Here the wave functions, time, and spatial coordi- 
nates are measured in the units of (h/muj)~ 3 / 4 , uj^ 1 , 
and (h/mtv) 1 / 2 , respectively, where u is the charac- 
teristic frequency of the trap for the atomic BEC 
fraction. The interaction strengths Uij are measured 
in the units of (huj)^ 1 (h/mijj)^ 3 ^ 2 , and the paramet- 
ric coupling strength, \, is measured in the units of 
{Tiuj)^ 1 (h/mus)~ 3 l A . We also introduce the dimension- 
less Raman detuning parameter 8, measured in the units 
of lo^ 1 , and the ratio between the trap frequencies for the 
molecular and atomic condensate fractions, Q — u) m /u). 

TABLE I. Model parameters for 87 Rb AMBEC. 



Model 


Parameter 


Dimensionless 


parameters 


values 


values 


U aa 


5.3 x lO^Jm 3 | 




0.062 


U am 


-6.9 x 10" 51 Jm 3 J 




-0.084(±0.07) 


Umm 


9.9 x 10~ 51 Jm 3 | 




0.12 


X 


7.4 x 10 _41 Jm 3/2 




1.09 



To specify the model parameters, we consider the case of 
87 Rb, for which the recent experiments on the formation 
of ultracold molecules were conducted . Realistic val- 
ues for the parameters can be deduced from a number of 
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available experimental data (see Table 1). We also take 
uj/2-k = 100, and £1 = 1. The atom-molecule two-body 
interaction strength was measured very recently p| . De- 
spite large uncertainty in the value of the corresponding 
s-wave scattering length, a am = — 180ao± 150ao, it is ap- 
parent that, unlike the atom- atom interaction, the atom- 
molecule interaction is attractive. This fact, as we show 
below, results in a competition of attractive and repulsive 
interactions in an AMBEC. 

To find possible steady states, we look for radially- 
symmetric solutions of Eqs. (|l|) in the form: i/) a (r,t) = 
(j) a (r) exp(-i/x a t), ip m (r, t) = (j> m {r) enp{-idt~i^ m t+i9), 
where 4>o,{t) and </> m (r) are positively defined, real func- 
tions, /x a and \i m are the chemical potentials, and 9 is 
a relative phase of the condensate fractions. The sub- 
stitution of the above ansatz into Eqs. (Q) reveals that, 
in order for a stationary AMBEC to exist, two condi- 
tions must be satisfied: \i m = 2fi a = 2/i and 9 — {0, 7r}. 
The first condition is reminiscent of the phase-matching 
of a fundamental and second-harmonics waves in an op- 
tical parametric process. The second condition implies 
that there are at least two physically different stationary 
configurations of the AMBEC, with the two components 
being "in-" or "out-of-phase" . The equations for (f> a {r) 
and (f> m (r) become: 

i (A r - r 2 ) cf) a - o-x4>a4>m + 

-(U aa 4>l + U am 4> 2 m )<i>a = 0, 

1 (2) 
- (A r - Ar 2 ) <p m + 5<j) m - cr|0 2 + 2n<j) m 

-(U m m4>m + Uam4> 2 a )4>m = 0, 

where A r = d 2 /dr 2 + (2/r)d/dr and a = cos 9 = ±1. In 
what follows, we assume near-resonant parametric cou- 
pling and set 6 = 0. 

From Eq. (f| it can be seen that the intra- and inter- 
species interactions act as repulsive or attractive poten- 
tials for the condensate components. Without the para- 
metric (Raman) coupling, the net effect of the cubic inter- 
actions is determined by the sign of AU = U aa U mm — U 2 m 
||. Since the atom-molecule interaction is attractive, 
the net-attractive and net-repulsive regimes occur when 
AU < and AU > 0, respectively. In addition, the 
quadratic Raman coupling is effectively repulsive for in- 
phase components (a = +1), and attractive for out-of- 
phase components (a = —1). Depending on the relative 
phase, 9, the system admits two types of stationary AM- 
BEC states. Typical families of stationary states in net- 
repulsive regime are shown in Fig. 1; we identify them by 
the mnemonics am + (for a = +1) or am_ (for a = — 1). 

The am + states are dominated by the repulsive cu- 
bic and parametric interactions, and exist only in the 
presence of a trap. Since the system (^]) admits station- 
ary solutions in the form of a pure molecular condensate, 



the family of the hybrid solutions branches off the pure 
molecular state at a bifurcation point (upper open circle 
in Fig. 1). In the vicinity of this point the atomic fraction 
is small. The physical meaning of this bifurcation and the 
coexistence of the pure molecular and AMBEC solution 
families is the following: unless some of the atomic con- 
densate is present, the two-component AMBEC will not 
form again, if a complete transformation into the molec- 
ular fraction has occured. This is a feature of the GP 
model which ignores spontaneous transitions from the 
molecular state (see discussions in p0|). 
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FIG. 1. Top: Families of AMBEC stationary states at 
Uam = —0.054: trapped out-of-phase (a2) and in-phase (b) 
states (solid), untrapped out-of-phase state (oi, dotted) and 
pure molecular BEC state (dashed). The trapped am l _ so- 
lution (solid) branches off the molecular BEC state at the 
bifurcation point (inset), (a) Radial wave functions, (j> a and 
4> m e' e , of trapped am^_ (dashed) and untrapped am"' (solid) 
states for fi = —1.45 and fi = —3.7, respectively, shown with 
the trapping potential V(r) = r 2 (dotted); N = 5.04 x 10 4 . 
(b) Radial wave functions of a trapped am+ stationary state 
for fi = 15 and N = 1.66 x 10 5 , shown with the Thomas-Fermi 
approximation for (f> m matching the number of atoms in the 
molecular fraction (dotted). 

Remarkably, the shape of the molecular component re- 
mains close to that obtained by the Thomas-Fermi ap- 
proximation even far from the bifurcation point, while 
the atomic component is moderate [see Fig. 1(b)]. This is 
due to the fact that, for the specific values of parameters 
corresponding to 87 Rb case (see Table 1), the effective po- 
tential experienced by the molecular component is dom- 
inated almost entirely by the molecular self-interaction; 
the inter-species and Raman interactions, being of dif- 
ferent sign and comparable magnitude, counteract and 
nearly cancel each other. Far from the bifurcation point, 
the density-dependent cubic interactions dominate over 
Raman coupling and the shape of the condensate deviates 
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from that found by the Thomas-Fermi approximation. 

In the presence of a trap, the am_ states branch off 
the pure molecular state at a small positive \i (Fig. 1, 
inset) and exist for /i < 0. For AC/ w 0, the cubic inter- 
actions nearly cancel each other so that the shapes of the 
condensate components are almost mirror images of each 
other, and the condensate is strongly affected only by 
the attractive Raman coupling. The resulting attractive 
interaction leads to self-trapping of the two-component 
condensate. Consequently, the am_ states can exist with 
(ami) or without (ami) the external trap. These states 
also exist when AU ^ 0, and the corresponding solutions 
for the net-repulsive regime are shown in Figs. l(a,b). 
Even a subtle change in the magnitude of AU results 
in a dramatic change of the peak density of the station- 
ary states. For a fixed AC/ and increasing particle num- 



ber, the 



solutions broaden spatially while keeping 



roughly the same peak density, asymptotically approach- 
ing the unbounded, homogeneous AMBEC. In contrast, 
the spread of the ami states is hindered by the exter- 
nal trap, so they rapidly grow in the peak density. As a 
consequence, the two initially close families of am_ de- 
part at large N (see Fig. 1). At the moderate values of 
N, while the families are still close, the condensate wave 
functions are weakly affected by the confining trap. 
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FIG. 2. Existence domains in the parameter space 
{ — U am , /J,}, (a) Trapped (hashed) and untrapped (dotted) 
stationary out-of-phase states, (b) Trapped in-phase AM- 
BEC (hashed) and trapped pure molecular BEC (dotted). 
Dash-dotted line marks the border AU = between the 
net-attractive (above the line, AU < 0) and net-repulsive 
(below the line, AU > 0) regimes. 

We note that the specific choice of parameters in Ta- 
ble 1 does not affect the generality of our results as the 
magnitude of AC/ only affects the densities of the station- 
ary AMBEC fractions. However, the properties of the 
AMBEC stationary states in the net-attractive regime 
of cubic interactions differ drastically from those in the 
net-repulsive regime. The two different regimes can be 
entered by varying either of the interaction strengths Uij 
(e.g., by tuning U aa through a Feshbach resonance). Here 
we vary the U am , within the large uncertainty in its value 
(see Table I). The existence regions for both in-phase 
and out-of-phase states are shown in Figs. 2(b) and 2(a), 
respectively. For the am + states, the prevailing attrac- 
tive interactions above the line AC/ = ramps up the 



density of the condensate at the center of trap, until the 
AMBEC states cease to exist at a certain C/ am (/z) [see 
Fig. 2(b)]. The ami states exist everywhere in the net- 
attractive domain [see Fig. 2(a)], and since both cubic 
and quadratic (parametric) interactions are effectively at- 
tractive, the AMBEC grows in central density, becoming 
highly localized with growing Both am + and ami 
states in a net-attractive regime are likely to be unstable 
against collapse. 

The AMBEC states in the net-repulsive regime are of 
great interest from the experimental point of view, pro- 
vided they are stable. A linear stability analysis for the 
homogeneous states, the amplidudes of which are asymp- 
totically approached by the ami* and ami states at the 
large particle number, show that both these states can 
be stable. Furthermore, our numerical analysis of the dy- 
namical stability of AMBEC states performed by solving 
the system (Q) with the initial conditions from different 
regions of the existence domain (see Fig. 2) confirmed 
that the am- states are likely to be stable, whereas all 
am+ states are dynamically unstable. 
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FIG. 3. (a) Evolution of an initially trapped, dynam- 
ically stable, stationary state ami, for /i = —2.5 and 
Uam = —0.054, released from the trap at t = 10 (dashed). 
(b,c) Evolution of the atomic and molecular components. Pa- 
rameters for the initial state are as in Table 1. 

Technically, the formation of an AMBEC should begin 
within a trap. However, the co-existence, in the same pa- 
rameter region, of two types of dynamically stable states, 
ami and am™', is a key to creation of AMBEC droplets 
without a trap. Namely, an ami state can be created, 
and then the trap can be turned off, which can cause a 
transition to the dynamically stable ami* state, provided 
the shapes and number of atoms of the condensates in the 
two states are close. This can be achieved by tuning the 
parameters so that the respective stationary states are 
weakly affected by the trapping potential (see Fig. 1). 
An example of the typical dynamics in this case is shown 
in Fig. 3, where the initially trapped out-of-phase state 
displays stable behavior. After the trap is turned off, the 
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AMBEC quickly settles into small-amplitude oscillations 
around a new (untrapped) stable stationary state. 

A natural route to the formation of a hybrid state 
is generation of the molecular component from a pure 
atomic condensate. Due to the parametric coupling, 
there is no stationary state corresponding to a pure 
atomic fraction, so formation of the molecular fraction 
starts immediately || . The initial stages of the AMBEC 
condensate formation are shown in Fig. 4(a). Without 
any dissipative processes, it is hard to anticipate that a 
dynamically stable stationary state would be created in 
a trap. However, even in this case the unstable dynamics 
of the condensate can be suppressed by turning the trap 
off. Figures 4(a-c) demonstrate that, after initial tran- 
sitional stage, the untrapped condensate does settle into 
oscillations around a stable state. The transition, how- 
ever, takes much longer than in the case of release of an 
already formed stationary am*! condensate (cf. Fig 3). 
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FIG. 4. (a) Generation of a trapped AMBEC from an 
atomic fraction (Gaussian initial profile) with subsequent re- 
lease from the trap at t — 2 (dashed line). Inset shows the 
long-term evolution. (b,c) Evolution of the atomic and molec- 
ular components. Parameters are as in Fig. 3 

It is important to mention that, in order to understand 
how a realistic AMBEC will relax to one of the stationary 
states, one should account for the various effects of losses. 
This can be done by adding phcnomcnological dissipa- 
tive terms into the coupled CP equations (see also Refs. 
HU). Additionally, one could take into account the effect 
of quantum fluctuations, by going beyond the two-mode 
approximation of the GP equations Jl|,|l^] , and deriving 
the stochastic field equations for the full quantum states 
of the condensates that include the GP equations as the 
noiseless limit [[To) . The resulting effective damping of the 
coherent atomic-molecular transfer oscillations will mod- 
ify the process of formation of the AMBEC (especially 
important after a near total conversion to the molecular 
fraction ]l0[]). However, all these effects will not signif- 
icantly affect the dynamically stable stationary states, 



once and if, the system has settled into one of them, and 
they can be easily accounted for in the framework of the 
adiabatic approximation. 

In conclusion, we have analyzed the mean-field model 
of an atomic-molecular condensate and identified the ex- 
istence domains for its stationary states. Our main find- 
ing is that, depending on the relative phase between the 
molecular and atomic fractions, two types of stationary 
states can be created in a confining trap. These are: (i) 
Thomas-Fermi-like states, existing due to a balance of 
the trapping potential and nonlinear interactions, and (ii) 
self-trapped states, dominated by the (attractive) para- 
metric inter-species interactions. In a certain region of 
the existence domain, the latter states are only weakly 
affected by the trap, and, upon release from the trap, 
they can be transformed into dynamically stable AM- 
BEC droplets that are localized even without a trap. We 
believe that this result is potentially useful for the atom 
(molecular) laser applications where the formation of sta- 
ble untrapped coherent matter waves is needed. 

The authors are indebted to P.D. Drummond, J. Hope, 
N. Robins, and C. Savage for useful discussions. 
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